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OPERATING A VEHICLE CONTROL SYSTEM 

5 

TECHNICAL FIELD OF THE INVENTION 

The present invention relates generally to motor vehicles. More particularly, the 
10 invention relates to a strategy for operating a vehicle stability control system. 

BACKGROUND OF THE INVENTION 

The dynamics of ground vehicles are largely determined by forces between its 
tires and the road. Knowledge of these forces is beneficial for most active vehicle control 

15 systems. Several strategies have been developed for estimating tire forces in three axes 
(e.g., longitudinal, lateral, and vertical). For example, it has been reported that the 
measurement of deformation of a magnetized tire may be used to determine estimates of 
tire longitudinal and lateral forces. This strategy may have several shortcomings 
including the need for transmission of data from a rotating tire, a requirement of an 

20 excessive amount of data processing, and a lack of robustness due to a need for using tire 
stiffness characteristics in relating tire forces to measured deformation. The tire stiffness 
characteristics may further be complicated by variations in pressure, rotation speed, and 
temperature. 

One strategy for estimating tire vertical forces using inertial sensors is disclosed 
25 in US patent 6,593,849 issued July 15, 2003 to Chubb, et al A two wheel lift off 

condition (i.e., when the vertical force on two tires of the same side of vehicle is zero) 
may be detected thereby indicating a danger of vehicle rollover. The strategy disclosed 
by the '849 patent, however, is limited to a calculation of vertical tire forces and 
therefore does not provide a determination of longitudinal and lateral forces. 

30 
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In US patent No. 6,549,842 issued April 15, 2003 to Hac, et al. 9 a strategy for 
estimating a surface coefficient of adhesion for each wheel is disclosed, which involves a 
determination of tire dynamic forces. This strategy, however, does not account for the 

5 effects of vehicle body pitch and roll dynamics or wheel vertical dynamics, which may 
limit the accuracy of estimates under dynamic conditions (e.g., driving on a rough road). 
Moreover, tire forces are calculated from models that depend on several variables (e.g., 
tire longitudinal slip, tire sideslip angle, and surface coefficient of adhesion), which 
cannot be measured directly (at a reasonable cost) and are difficult to estimate with 

10 precision. 

Therefore, it would be desirable to provide a strategy for operating a vehicle 
control system based on determined dynamic tire forces that overcomes the 
aforementioned and other disadvantages. 

SUMMARY OF THE INVENTION 

A first aspect of the present invention provides a method of determining dynamic 
tire forces. A plurality of measured vehicle variables is determined. At least one 
suspension force is determined based on at least one of the determined measured vehicle 
variables. At least one tire vertical force is determined based on the at least one 
determined suspension force and at least one of the determined measured vehicle 
variables. At least one tire longitudinal force is determined based on at least one of the 
determined measured vehicle variables. At least one axle lateral force is determined 
based on at least one of the determined measured vehicle variables. At least one tire 
lateral force is determined based on the at least one determined axle lateral force and the 
at least one determined tire vertical force. At least one control component of the vehicle 
control system is modulated based on the at least one determined tire vertical force, the at 
least one determined tire longitudinal force, and the at least one determined tire lateral 
force. 
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A second aspect of the invention provides a computer usable medium including a 
program for determining dynamic tire forces. The computer usable medium includes 
computer readable program code for: determining a plurality of measured vehicle 

5 variables; determining at least one suspension force based on at least one of the 

determined measured vehicle variables; determining at least one tire vertical force based 
on the at least one determined suspension force and at least one of the determined 
measured vehicle variables; and determining at least one tire longitudinal force based on 
at least one of the determined measured vehicle variables. The computer usable medium 

10 further includes computer readable program code for: determining at least one axle lateral 
force based on at least one of the determined measured vehicle variables; determining at 
least one tire lateral force based on the at least one determined axle lateral force and the 
at least one determined tire vertical force; and modulating at least one control component 
of the vehicle control system based on the at least one determined tire vertical force, the 

15 at least one determined tire longitudinal force, and the at least one determined tire lateral 
force. 

A third aspect of the invention provides a system for determining dynamic tire 
forces. The system includes means for: determining a plurality of measured vehicle 
variables; determining at least one suspension force based on at least one of the 

20 determined measured vehicle variables; determining at least one tire vertical force based 
on the at least one determined suspension force and at least one of the determined 
measured vehicle variables; and determining at least one tire longitudinal force based on 
at least one of the determined measured vehicle variables. The system further includes 
means for: determining at least one axle lateral force based on at least one of the 

25 determined measured vehicle variables; determining at least one tire lateral force based 
on the at least one determined axle lateral force and the at least one determined tire 
vertical force; and modulating at least one control component of the vehicle control 
system based on the at least one determined tire vertical force, the at least one determined 
tire longitudinal force, and the at least one determined tire lateral force. 
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The foregoing and other features and advantages of the invention will become 
further apparent from the following detailed description of the presently preferred 
embodiments, read in conjunction with the accompanying drawings. The detailed 
5 description and drawings are merely illustrative of the invention, rather than limiting the 
scope of the invention being defined by the appended claims and equivalents thereof. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic view of a vehicle including a system for operating a vehicle 
10 control system in accordance with the present invention; 

FIG. 2 is a flow diagram of an algorithm for determining dynamic tire forces in 
accordance with the present invention; 

FIG. 3A illustrates a quarter-car model with forces acting in a vertical direction 
on a vehicle body and wheel; 
15 FIG. 3B illustrates lateral forces acting on a vehicle represented as a bicycle 

model; 

FIG. 4 illustrates forces and torques affecting wheel rotation of a vehicle; 
FIG. 5 illustrates a simplified diagram of forces acting upon a vehicle body and 
one wheel; 

20 FIG. 6 illustrates a detailed full vehicle model comparison of simulated tire 

forces to estimations derived from the algorithm of FIG. 2; and 

FIG. 7 illustrates a comparison of simulated vertical tire force to estimations 
derived from the algorithm of FIG. 2 during driving on a rough road surface. 
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DESCRIPTION OF THE PRESENTLY PREFERRED EMBODIMENTS 

Referring to the drawings, wherein like reference numerals refer to like elements, 
FIG. 1 is a schematic view of a vehicle, shown generally by numeral 10, with a system 
5 20 for determining dynamic tire forces in accordance with the present invention. Those 
skilled in the art will recognize that the vehicle 10 and system 20 may include a number 
of alternative designs and may be employed in a variety of applications. For example, as 
will be described below, the vehicle 10 may include various sensor "packages" as part of 
different system 20 embodiments. 

10 In one embodiment, the vehicle 10 may include two axles 30, 32, each flanked by 

wheels 34, 36, 38, 40. Vehicle 10 may further include a controller 50 including a digital 
microprocessor 52 for executing an algorithm for determining dynamic tire forces in 
accordance with the present invention. Controller 50 may be programmed to process a 
plurality of input signals received from a sensor array 60 and a vehicle control system 

1 5 (VCS) 70. Further, the controller 50 may provide output signals to the sensor array 60 
and VCS 70, which may be used for modulating (i.e., enhance or otherwise control in an 
advantageous manner) the performance of various chassis control algorithms. In one 
embodiment, the VCS 70 may modulate various systems including, but not limited to, a 
chassis controllable suspension, an anti-lock braking system (ABS), a traction control 

20 system (TCS), a vehicle stability enhancement system (VSE), an active rear steering 

system, and an anti-rollover system. Those skilled in the art will appreciate that the VCS 
70 may modulate numerous factors influencing vehicle comfort, handling, stability, and 
various other factors based on the determined tire forces. 

The methods, algorithms, integration of data, and determinations (e.g., 

25 calculations and estimations), of the presently preferred embodiments, including those 
based on equations or value tables, may be performed by a device such as the 
microprocessor 52. Additional sensors 42, 44, 46, 48 may be optionally positioned 
within or associated with the wheels 34, 36, 38, 40 and/or vehicle 10 suspension for 
determining wheel 34, 36, 38, 40 speed, vertical acceleration, and other information. 

30 This information may be communicated to the sensor array 60 as denoted by the arrows 
82, 84, 86, 88. 
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The computer usable medium including a program for operating the VCS 70 and 
program code associated with the presently preferred embodiments may be read into and 
stored in a memory portion 54 (e.g., ROM, RAM, and the like) for access by the 
5 microprocessor 52, as understood in the art. Furthermore, value tables, variables, 

parameters, and other numerical data may be stored as necessary in the memory portion 
54. Analog signal processing may be provided for some of the input signals received by 
the controller 50. For example, the information sent from the sensor array 60 may be 
low-pass filtered through analog filter(s) to reduce signal noise. 

10 For simplicity, the various sensors comprising the sensor array 60 are shown as a 

composite unit positioned at or near a center of gravity of the vehicle 10. Depending on 
the nature of the vehicle 10, the sensor array 60 may comprise a variety of individual 
sensors including, but not limited to, those sensing longitudinal, lateral, and vertical 
(a.k.a. normal) body accelerations, steering angle, roll rate, pitch rate, yaw rate, and 

15 wheel brake and drive torques (in addition to the wheel speed and vertical wheel 

acceleration information received by the wheel sensors 42, 44, 46, 48). Those skilled in 
the art will recognize that the various sensors may be substituted, added, or reconfigured 
to determine dynamic tire forces. The inventor contemplates a wide variety of sensor 
"packages" that may be used to effectively determine dynamic tire force in accordance 

20 with the present invention. Further, the positioning of the individual sensors and/or 
sensor array 60 may vary within the vehicle 10 and are not necessarily limited to 
placement at or near the vehicle 10 center of gravity as shown. 

FIG. 2 is a flow diagram of an algorithm for determining dynamic tire forces in 
accordance with the present invention. The algorithm may be implemented by the 

25 controller 50 shown in FIG. 1 and may provide means for estimating tire forces in 

longitudinal, lateral, and vertical directions using primarily inertial sensors located in the 
sensor array 60 and (optionally) the wheels 34, 36, 38, 40. In addition, wheel speed and 
estimates of brake and drive torques may be used primarily in determination of tire 
longitudinal forces. The algorithm may be derived from equations of vehicle 10 vertical 

30 and lateral dynamics and the wheel 34, 36, 38, 40 rotary dynamics. In addition to 
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measured variables, vehicle 10 parameters may be used. As some of the inertial 
parameters of the vehicle 10 may vary to some degree with payload, additional parameter 
adaptation algorithms may be implemented to improve accuracy of the algorithm. 
5 Tire vertical forces may be determined using measured longitudinal, lateral, and 

vertical accelerations of the vehicle 10 body, and the roll, pitch, and yaw rates. In 
addition, an estimate of vehicle 10 speed may be used to determine the influence of 
aerodynamic lift force and pitch moment. If a higher accuracy is desired, vertical 
accelerations of the wheels 34, 36, 38, 40 or relative velocities/positions of the wheels 34, 

10 36, 38, 40 with respect to the vehicle 10 body may also be used, depending on which 

sensors are available. Longitudinal tire forces may be determined from wheel 34, 36, 38, 
40 speeds, estimates of brake and driving torques, and, optionally, tire vertical forces. 
The tire lateral forces may be determined using primarily lateral acceleration, yaw rate, 
and tire vertical forces, and, optionally, tire longitudinal forces and steering angle, if the 

1 5 latter are available. 

In one embodiment, the algorithm provides a strategy for estimating dynamic tire 
forces in three axes (e.g., longitudinal, lateral, and vertical). Those skilled in the art will 
recognize that the aforementioned three directions are relative and may vary slightly due 
to the orientation of the vehicle 10, road, wheels 34, 36, 38, 40, etc. Ideally, the 

20 estimation of the dynamic tire forces may be obtained from relatively inexpensive 

sensors as part of the sensor array 60 and, optionally, associated with wheel 34, 36, 38, 
40 drive-train and brake systems. This information may be obtained from inertial sensors 
(e.g., accelerometers and rate sensors), and in the case of brake torque, their estimates 
may be provided by brake and power train controllers (not shown). The algorithm is at 

25 least partially based on the knowledge of vehicle 10 inertial parameters (e.g., masses and 
moment of inertia) and the location of the vehicle 10 center of gravity. 
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For most vehicles, these parameters do not vary much with payload. However, 
the accuracy of the algorithm may be improved if the vehicle payload condition is 
estimated via a parameter adaptation algorithm. Such adaptation algorithms have been 
5 previously disclosed. For example, in US patent application 10/440,604 filed May 19, 
2003 by Hac et ai, a recursive algorithm is disclosed for estimating vehicle payload 
condition using information from roll rate and lateral acceleration sensors (or pitch rate 
and longitudinal acceleration information). In US patent 6,424,907 issued July 23, 2003 
to Rieth, et al, corrections to a vehicle center of gravity location (from the nominal 

10 values) are determined using information from longitudinal and lateral acceleration and 
static suspension deflection sensors or pressure sensors in air spring of self-leveling 
suspensions (if available). 

Several key concepts of the algorithm are illustrated in FIGS. 3A and 3B. FIG. 
3A illustrates a quarter-car model with forces acting in a vertical direction on a vehicle 

1 5 body and wheel. Here m hc represents the body (spaing) mass (per corner), m w the un- 
sprung (wheel) mass, F s is the suspension force, F z is the tire vertical force and F m 
represents the forces acting directly on the vehicle body (other than the suspension force), 
such as inertial forces due to braking/accelerating or cornering and aerodynamic lift 
force. The equations of motion of the body and the wheel may be expressed as: 

20 



25 



where a zb and a 2W are the measured accelerations of the body and the wheel, respectively. 
They include components of gravity acceleration. That is, on level road: 



<** = ' z b-g (2 ) 
a = z — 2 
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where z A and z w are the vertical accelerations of body and wheel, respectively. The 
equations (1) can be solved for an unknown vertical load, F 2 . This yields: 

5 F z = F m -m hc a zh -m w a :w (3) 

When the vehicle is at rest on level road, F m = 0 and a z b - a zw = - g. 
Consequently, F z ={ir\ } +n\)g, which is the vertical force, is equal to the static load. Thus, 

the vertical force calculated from equation (3) includes the static load, as well as dynamic 
10 loads due to inertial forces resulting from road inputs and vehicle maneuvers. Also, the 
effect of road inclinations on vertical force is included in measured vertical accelerations. 
There is no requirement to measure suspension force. This concept may be extrapolated 
to a full car model. 

If an accelerometer measuring wheel vertical acceleration is not available, the last 
15 term in equation (3) may be ignored and the quarter body mass replaced by a quarter 
vehicle mass, to account for the weight of the wheel. This yields: 

F z = F m - m vc a zh (3b) 

20 If, however, a measurement of relative velocity between the wheel and the body is 

available (this can be measured directly or obtained by differentiation of suspension 
deflection measurement), v rc i = v zw - v zb , then the wheel acceleration may be expressed 
as: 

25 a zw =v rel +a zh (4) 

and this value substituted into equation (3). 
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FIG. 3B illustrates lateral forces acting on a vehicle represented as a bicycle 
model. Assuming a small steering angle, the equations of motion in the yaw plane may 
be expressed as: 



ma y = F yf + F yr 

l J>z = F yf a ~ F yr b 



where m is the total mass of vehicle, I zz is vehicle moment of inertia about the vertical 
axis, F v f and F yr are the front and rear axle lateral forces per axle, a y is measured lateral 
10 acceleration (which includes the gravity component due to road inclination), co z is the 
yaw acceleration, and a and b are the distances of vehicle center of mass to front and rear 
axle, respectively. The equations (5) may be solved for the lateral forces (per axle), 
yielding: 



mba v + I„(o„ maa v - I^co y 
F = 22 2 F = - (6) 



where L = a+b is vehicle wheelbase. On a uniform surface, the lateral forces may be 
distributed between left and right tires approximately in proportion to the corresponding 
vertical forces. If a higher accuracy is desired, the effects of differential braking and 
20 steering angle may be incorporated into equations (5) and (6). This may require 
measurement of the steering angle and tire longitudinal forces. 

The forces and torques affecting wheel rotation are illustrated in FIG. 4. The 
equation of motion may be expressed as: 

25 I W G) W = T d -T h - eF z - F x r w (7) 
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where / w is a wheel rotary moment of inertia, a** is the wheel rotational speed, T c i is the 
drive torque applied to the wheel, T b is the brake torque, F z is the vertical force, F x is the 
longitudinal force, r w is the tire radius, and e is the offset of vertical force from the wheel 
5 centerline. If the drive and/or brake torques are estimated (e.g., by power train or brake 
controllers) and the wheel acceleration is estimated from wheel speed, the tire 
longitudinal force may be determined from equation (7) as a function of known variables: 



F = T d -T b -eF 2 -I w d) w (g) 



10 



It should be noted that the determination of lateral forces per tire and longitudinal 
forces may depend on the vertical tire forces, although in the latter case this dependency 
is weak, as the term eF z is usually small. However, calculation of the longitudinal forces 
does not depend on the lateral forces and calculation of lateral forces may be only weakly 
15 dependent on longitudinal forces (only in the case when differential braking is taken into 
account). Thus, the order of calculating longitudinal and lateral tire forces may be 
arbitrary. 

In the algorithm, the determination of dynamic tire lateral forces may be 
extrapolated to the full car model. FIG. 5 illustrates a simplified diagram of forces 

20 acting upon a vehicle body and one wheel. The vehicle body may be subjected to four 
(or more or less) suspension forces and a number of forces and moments (comprising 
aerodynamic and inertial forces and moments), which in the quarter car model were 
lumped as a single force, F m . Specifically, the following forces may act on vehicle body: 
four suspension forces in four corners; gravity force; aerodynamic lift force, Fun; 

25 aerodynamic pitch moment, M p ; longitudinal, lateral, and vertical inertial forces, ma x , 
ma y , ma z ; and roll, pitch and yaw inertial moments, l xx d) x , ^ yy co y , I 2Z <*V All these forces 

and moments are now included in the described algorithm. 

Referring again to FIG* 2, in one embodiment, the algorithm may be divided into 
the following steps: 

30 
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A plurality of measured vehicle variables is determined (step 101). The variables 
may comprise receiving and filtering the following measured sensor information: body 
longitudinal acceleration, a x , body lateral acceleration, a y , body vertical acceleration, a Z b> 
5 wheel vertical accelerations, a ZW LF> a ZW RF, a 2W LR, a ZW RR, body roll rate, co x , body pitch rate, 
coy, and yaw rate, oo x , wheel speeds cd w lf> co wL r, co w rf, co w rr. The roll, pitch, and yaw 
rates may be differentiated (i.e., in practice, high-pass filtered) to obtain the roll, pitch, 
and yaw accelerations. Also, the wheel speeds may be differentiated to obtain wheel 
angular accelerations. 

10 It may be assumed that the body accelerations are determined by accelerometers 

placed at or close to a vehicle center of gravity. If accelerometers measuring wheel 
vertical accelerations are not available, but wheel to body relative velocities are obtained 
from suspension sensors, then the wheel vertical accelerations may be determined by 
adding the time derivatives of wheel relative velocities to the vertical accelerations of 

1 5 body corners at the wheel location. The latter values may be obtained from the measured 
vertical acceleration at the body center of mass and the roll and pitch rates. This may 
yield the following expressions for wheel vertical accelerations: 

<W = a :h ~ + K®x 1 2 + W (9a) 
20 a 2wRF = a zh - a h 6) y - t w 6) x 1 2 + v relRF (9b) 

= a zt> + b b®y + / 2 + KeUR ( 9c ) 

***R = a :b + b h d) y - t w 6) x 1 2 + v rem (9d) 

where a h and b b are the distances of vehicle body center of gravity to the front and rear 
25 axles, respectively, and t w is the track width (assumed the same for front and rear axles). 
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At least one suspension force is determined based on at least one of the measured 
vehicle variables (step 102). Suspension forces may be determined using linear and 
rotational vehicle body accelerations and optionally the estimates of aerodynamic lift 
5 force and pitch moment. The suspension forces (in the vertical direction) may be 
calculated as follows: 



10 
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20 where m b is the total sprung (body) mass, and Iyy, are the body moments of inertia 
about the roll and pitch axes, respectively, h is the height of vehicle body center of 
gravity above ground, h r is the height of body center of gravity above the roll axis, a& and 
bh are distances of vehicle body center of gravity from front and rear axle, respectively, L 
= ai, + bb is the vehicle wheel-base, and t w is the track width (it is assumed the same at 

25 front and rear). Furthermore, tj represents the ratio of rear to front suspension roll 
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stiffness (including suspension springs and roll bars). M p and F/,y, are the estimated 
aerodynamic pitch moment and lift force, respectively, which may be determined in a 
manner known to those skilled in art using vehicle speed and known aerodynamic pitch 
5 and lift coefficients for a particular vehicle. The variables a x , a y , and a z b denote the 

measured body accelerations in the longitudinal, lateral, and vertical directions, including 
the effects of gravity components. The variables co x , co y , and co z are the rates of rotation 
of vehicle body about the roll, pitch, and yaw axes, wherein the dot denotes a derivative 
with respect to time. In practice, time derivatives may be approximated by high pass 

10 filters. Note that the effects of road inclinations on suspension forces may be accounted 
for in measured accelerations. 

The vehicle mass and moments of inertia used in the above equations may 
represent the nominal values; these may correspond, for example, to a vehicle with a 
driver and one passenger. However, these parameters may optionally be updated by 

1 5 determining a parameter adaptation that estimates vehicle loading conditions (step 102b). 

At least one tire vertical force is determined based on the suspension force and at 
least one of the measured vehicle variables (step 103). The vertical forces may be 
determined using the suspension forces determined in step 102 and measured vertical 
accelerations of the wheels expressed as follows: 

20 



FzLF = ~ F sLF 


(11a) 


F zRF ~ ~ F sRF ~ m w f a zwRF 


(Ub) 


F zLR ~ ~ F sLR ~~ m wr a zwLR 


(11c) 


F-RR = ~FsRR ~ m wr a zwRR 


(lid) 



25 

where m w f and m nv denote the unsprung masses of front and rear wheel assemblies, and 
a zw jj refer to the measured vertical accelerations of LF, RF, LR and RR wheels. These 
accelerations include the gravity acceleration. 
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If wheel accelerations are not available, then the last terms in equations (11) may 
be deleted and the vertical forces are equal to negatives of suspension forces (i.e., the 
negative sign is due to assumed sign convention). However, in this case, the body mass 
5 m h in equations (10) may be replaced by the total vehicle mass, m, and the distances h, h r , 
a, and b defined with respect to vehicle center of gravity, rather than the vehicle body 
center of gravity (in practice, the difference in location of both centers is small). 

At least one tire longitudinal force is determined based on at least one of the 
measured vehicle variables (step 104). Longitudinal tire forces may be determined using 
10 measured wheel speeds and estimated driving, braking, tire rolling resistance torques. 
The equations used for this purpose may be expressed as follows: 

xLF = \ 1 l *> 

„ T dRF - T hRF - eF zRF - I W fd) wRF /nM 

xRF = ( 1 2b ) 

15 f = ^ dLR ~^*>LR ~ e FzLR ~ ^wr^wLR (12c) 

xLR r 

wr 

f - ^ ilRR ~ ^ hRR ~~ e ^z RR ~ ^wr^wRR ( 1 2d) 

XRR y 

wr 



where T d is the estimated drive torque, T b denotes the estimated brake torque, eF z 
represents the moment due to rolling resistance, I w is the rotary inertia of the wheel (f- 
20 front, r-rear), ^ is the measured wheel rotational speed, and r w is the wheel radius. Note 
that the rolling resistance moment is a function of vertical load, but may also depend on 
vehicle speed, or may be a constant if lesser accuracy is acceptable. 
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At least one axle lateral force is determined based on at least one of the measured 
vehicle variables (step 105). Lateral forces per axle may be determined using the 
measured lateral acceleration and roll rate and known vehicle parameters according to: 



mba v + I. z cb. 
~L 



F y/ = (13a) 



F= L—^ (13b) 



where m is the total mass of vehicle, I zz is vehicle moment of inertia about the vertical 
1 0 axis, F y f and F yr are the front and rear axle lateral forces, a y is measured lateral 

acceleration (which includes the gravity component due to road inclination),^ is the 
yaw acceleration (obtained through differentiation of measured yaw rate), a and b are the 
distances of vehicle center of mass to front and rear axle, respectively, and L = a + b is 
vehicle wheelbase. The accuracy of the above equations may be improved by including 
1 5 the effect of differential braking and the effect of steering angle if the steering angle 
measurement is available. In this case the equations may be expressed as follows: 

p = mba y +L z <b z -(F xU r+F xRF y>sm8-^ 
yf " ~~ Lcos£ 

20 (14a) 

maa v - 1 zz 6) z - (F xLF + F xRF )a sin 8 + 0.5[(F xLF - F xRF )cosS + {F xLR - F xRR 



F.„ = 



LcosS 



(14b) 



25 
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At least one tire lateral force is determined based on the axle lateral force and the 
tire vertical force (step 106). Individual tire lateral forces may be determined using the 
lateral forces per axle determined in step 105 and the vertical tire forces determined in 
5 step 103. 

F F FA F :of\A F :of- F :L F )\ ^ 

yLF ~ * F UF [F zof +e(F sof -F !LF )\ + F !KF [F !of + s(F !of — F zRF jj 



F. RF \F_„ f + e(F, - F., F )| 
10 F vRF = F„ ^-rr / „ ~ u „ i„ /„ (15b) 



F vLR = F -—^ ,„ 1 „ * V ? r r J L * n f „ — ( 1 5c) 



F; LF [ 


F !of +e(F zo/ - 


F !LF )\+F :HF [F :o/+ e(F :o/ 






F !LR [ 


F or +s{F :or -F !LR )] 




f :L A> 


*; 0r +e{F !or - 


F :LR )]+ F :HR [F zor + s(F :or 


-F :RR )] 




F : rr [ 


F zor +e{F :or -F zRR )] 




F :L «{> 


r +e(F - 

zor \ zor 


F :LR )]+ F :RR [F :or + s{F !or 


-F :RR )} 



F vRR = F w rr-rr: , J RRl fr vT l° r r- — ( 1 5d) 



15 

where F zo /and F zor are the nominal static tire vertical forces for a front and rear tires, 
respectively. The coefficient, 8, is a constant that describes sensitivity of tire lateral 
forces to the changes in vertical load. A typical value of e may be about 0.3. 

Note that when £ = 0, the tire lateral forces on the left and right side of vehicle 
20 become proportional to vertical loads. That is, equations (15) simplify to: 
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F + F 

zLF ^ r zRF 



(16b) 



F 9 



zLR 



^zLR + ^zRR 



FyRR ~ F> 



■RR 



F + F 

1 zLR T ' z/?/? 



(16c) 



(16d) 



The equations (16) are the simplified versions of equations (15) and may be used when 
less accuracy is acceptable or when the tire sensitivity coefficient, 8, is small. 

10 At least one control component of the VCS is modulated based on the tire forces 

determinations (e.g., the tire vertical, longitudinal, and lateral forces; step 107). In one 
embodiment, the VCS may modulate one or more control component algorithm(s) of a 
chassis controllable suspension, an ABS, a TCS, a VSE, an active rear steering system, 
and an anti-rollover system based on the determined tire force determinations. In another 

1 5 embodiment, the VCS may modulate one or more control component algorithm(s) 
influencing vehicle comfort, handling, stability, or various other factors based on the 
determined tire forces. One skilled in the art will understand modulation of one or more 
control components of a VCS. 

The algorithm described herein relies on the knowledge of vehicle parameters 

20 related to inertial properties and the location of vehicle center of mass. These parameters 
may vary to some extent with vehicle payload. The robustness of the algorithm described 
herein may further be improved if the parameters such as mass, moments of inertia, and 
location of center of mass are adapted to vehicle payload. 
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An example of a method that permits estimation of vehicle payload using roll 
rate, pitch rate, lateral, and longitudinal acceleration sensors and may be adapted for use 
with the present invention is disclosed in United States Patent Application 10/440,604, 
5 filed May 19, 2003, the entirety of which is hereby incorporated by reference. In this 
method, a payload parameter, w, is determined. The value of this vehicle payload 
parameter varies from 0 to 1 , with 0 corresponding to a vehicle with a driver only and 1 
corresponding to a vehicle with the driver and a full payload. Vehicle parameters used 
herein, such as mass, moments of inertia and distances of center of gravity, may be either 

10 known or may be determined for the two extreme cases of loading, that is, for the vehicle 
with the driver only and for the vehicle with the driver and full payload. Consequently, 
each one of the vehicle actual parameters may be updated during vehicle operation by 
using the parameter, w, to interpolate between the two extreme values. For example, if 
vehicle body mass, m b , is nibempty for a vehicle with driver only and m b ioaded for vehicle 

1 5 with full payload, then the actual body mass may be determined as follows: 



»h = "W (1 " w) + m uoade<t w ( 1 7 ) 

Analogous equations apply to the other parameters. 

20 FIG. 6 illustrates a detailed full vehicle model comparison of simulated (solid 

lines) tire forces to estimations (dashed lines) derived from the algorithm. The results 
were obtained for a cornering and braking maneuver performed on a smooth, dry road 
surface. Steering was initiated half a second into the maneuver and braking at three 
seconds. The resulting tire forces in three directions for left front and right rear tires are 

25 shown. 

FIG. 7 illustrates a comparison of simulated (solid lines) vertical tire force to 
estimations (dashed lines) derived from the algorithm during driving on a rough road 
surface. A sinusoidal road input of frequency 3 Hertz and amplitude of 2.5 cm (about 1 
inch) was applied to the left wheel. . 

30 
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While the embodiments of the invention disclosed herein are presently considered 
to be preferred, various changes and modifications may be made without departing from 
the spirit and scope of the invention. For example, the vehicle and method of 

5 determining dynamic tire forces are not limited to any particular design, configuration, 
sequence, or arrangement. The sensors, and controller configuration, size, shape, 
geometry, location, orientation, number, and function may vary without limiting the 
utility of the invention. Furthermore, the method may be accomplished by numerous 
alternative strategies, may include additional steps, and vary in step order. 

10 Upon reading the specification and reviewing the drawings hereof, it will become 

immediately obvious to those skilled in the art that myriad other embodiments of the 
present invention are possible, and that such embodiments are contemplated and fall 
within the scope of the presently claimed invention. The scope of the invention is 
indicated in the appended claims, and all changes that come within the meaning and 

1 5 range of equivalents are intended to be embraced therein. 
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